Search for: All records

ABSTRACT The initial mass and metallicity of stars both have a strong impact on their fate. Stellar axial rotation also has a strong impact on the structure and evolution of massive stars. In this study, we exploit the large grid of GENEC models, covering initial masses from 9 to 500 $${\rm M}_{\odot }$$ and metallicities ranging from $$Z=10^{-5}$$ (nearly zero) to 0.02 (supersolar), to determine the impact of rotation on their fate across cosmic times. Using the carbon–oxygen core mass and envelope composition as indicators of their fate, we predict stellar remnants, supernova engines, and spectroscopic supernova types for both rotating and non-rotating stars. We derive rates of the different supernova and remnant types considering two initial mass functions to help solve puzzles such as the absence of observed pair-instability supernovae. We find that rotation significantly alters the remnant type and supernova engine, with rotating stars favouring black hole formation at lower initial masses than their non-rotating counterparts. Additionally, we confirm the expected strong metallicity dependence of the fates with a maximum black hole mass predicted to be below 50 $${\rm M}_{\odot }$$ at SMC or higher metallicities. A pair-instability mass gap is predicted between about 90 and 150 $${\rm M}_{\odot }$$, with the most massive black holes below the gap found at the lowest metallicities. Considering the fate of massive single stars has far-reaching consequences across many different fields within astrophysics, and understanding the impact of rotation and metallicity will improve our understanding of how massive stars end their lives, and their impact on the Universe. 

ABSTRACT One-dimensional (1D) stellar evolution models are widely used across various astrophysical fields, however they are still dominated by important uncertainties that deeply affect their predictive power. Among those, the merging of independent convective regions is a poorly understood phenomenon predicted by some 1D models but whose occurrence and impact in real stars remain very uncertain. Being an intrinsically multi-D phenomenon, it is challenging to predict the exact behaviour of shell mergers with 1D models. In this work, we conduct a detailed investigation of a multiple shell merging event in a 20 M$$_\odot$$ star using 3D hydrodynamic simulations. Making use of the active tracers for composition and the nuclear network included in the 3D model, we study the merging not only from a dynamical standpoint but also considering its nucleosynthesis and energy generation. Our simulations confirm the occurrence of the merging also in 3D, but reveal significant differences from the 1D case. Specifically, we identify entrainment and the erosion of stable regions as the main mechanisms that drive the merging, we predict much faster convective velocities compared to the mixing-length theory velocities, and observe multiple burning phases within the same merged shell, with important effects for the chemical composition of the star, which presents a strongly asymmetric (dipolar) distribution. We expect that these differences will have important effects on the final structure of massive stars and thus their final collapse dynamics and possible supernova explosion, subsequently affecting the resulting nucleosynthesis and remnant. 

ABSTRACT The treatment of convection remains a major weakness in the modelling of stellar evolution with one-dimensional (1D) codes. The ever-increasing computing power makes now possible to simulate in three-dimensional (3D) part of a star for a fraction of its life, allowing us to study the full complexity of convective zones with hydrodynamics codes. Here, we performed state-of-the-art hydrodynamics simulations of turbulence in a neon-burning convective zone, during the late stage of the life of a massive star. We produced a set of simulations varying the resolution of the computing domain (from 1283 to 10243 cells) and the efficiency of the nuclear reactions (by boosting the energy generation rate from nominal to a factor of 1000). We analysed our results by the mean of Fourier transform of the velocity field, and mean-field decomposition of the various transport equations. Our results are in line with previous studies, showing that the behaviour of the bulk of the convective zone is already well captured at a relatively low resolution (2563), while the details of the convective boundaries require higher resolutions. The different boosting factors used show how various quantities (velocity, buoyancy, abundances, and abundance variances) depend on the energy generation rate. We found that for low boosting factors, convective zones are well mixed, validating the approach usually used in 1D stellar evolution codes. However, when nuclear burning and turbulent transport occur on the same time-scale, a more sophisticated treatment would be needed. This is typically the case when shell mergers occur. 

ABSTRACT Our knowledge of stellar evolution is driven by one-dimensional (1D) simulations. 1D models, however, are severely limited by uncertainties on the exact behaviour of many multidimensional phenomena occurring inside stars, affecting their structure and evolution. Recent advances in computing resources have allowed small sections of a star to be reproduced with multi-D hydrodynamic models, with an unprecedented degree of detail and realism. In this work, we present a set of 3D simulations of a convective neon-burning shell in a 20 M⊙ star run for the first time continuously from its early development through to complete fuel exhaustion, using unaltered input conditions from a 321D-guided 1D stellar model. These simulations help answer some open questions in stellar physics. In particular, they show that convective regions do not grow indefinitely due to entrainment of fresh material, but fuel consumption prevails over entrainment, so when fuel is exhausted convection also starts decaying. Our results show convergence between the multi-D simulations and the new 321D-guided 1D model, concerning the amount of convective boundary mixing to include in stellar models. The size of the convective zones in a star strongly affects its structure and evolution; thus, revising their modelling in 1D will have important implications for the life and fate of stars. This will thus affect theoretical predictions related to nucleosynthesis, supernova explosions, and compact remnants. 

ABSTRACT Our understanding of stellar structure and evolution coming from one-dimensional (1D) stellar models is limited by uncertainties related to multidimensional processes taking place in stellar interiors. 1D models, however, can now be tested and improved with the help of detailed three-dimensional (3D) hydrodynamics models, which can reproduce complex multidimensional processes over short time-scales, thanks to the recent advances in computing resources. Among these processes, turbulent entrainment leading to mixing across convective boundaries is one of the least understood and most impactful. Here, we present the results from a set of hydrodynamics simulations of the neon-burning shell in a massive star, and interpret them in the framework of the turbulent entrainment law from geophysics. Our simulations differ from previous studies in their unprecedented degree of realism in reproducing the stellar environment. Importantly, the strong entrainment found in the simulations highlights the major flaws of the current implementation of convective boundary mixing in 1D stellar models. This study therefore calls for major revisions of how convective boundaries are modelled in 1D, and in particular the implementation of entrainment in these models. This will have important implications for supernova theory, nucleosynthesis, neutron stars, and black holes physics. 

Abstract We present the state-of-the-art single-zone nuclear reaction networkWinNet, which is capable of calculating the nucleosynthetic yields of a large variety of astrophysical environments and conditions. This ranges from the calculation of the primordial nucleosynthesis, where only a few nuclei are considered, to the ejecta of neutron star mergers with several thousands of involved nuclei. Here we describe the underlying physics and implementation details of the reaction network. We additionally present the numerical implementation of two different integration methods, the implicit Euler method and Gears method, along with their advantages and disadvantages. We furthermore describe basic example cases of thermodynamic conditions that we provide together with the network and demonstrate the reliability of the code by using simple test cases. With this publication,WinNetwill be publicly available and open source at GitHub and Zenodo. 

Context. A realistic parametrization of convection and convective boundary mixing in conventional stellar evolution codes is still the subject of ongoing research. To improve the current situation, multidimensional hydrodynamic simulations are used to study convection in stellar interiors. Such simulations are numerically challenging, especially for flows at low Mach numbers which are typical for convection during early evolutionary stages. Aims. We explore the benefits of using a low-Mach hydrodynamic flux solver and demonstrate its usability for simulations in the astrophysical context. Simulations of convection for a realistic stellar profile are analyzed regarding the properties of convective boundary mixing. Methods. The time-implicit Seven-League Hydro (SLH) code was used to perform multidimensional simulations of convective helium shell burning based on a 25  M ⊙ star model. The results obtained with the low-Mach AUSM + -up solver were compared to results when using its non low-Mach variant AUSM B + -up. We applied well-balancing of the gravitational source term to maintain the initial hydrostatic background stratification. The computational grids have resolutions ranging from 180 × 90 2 to 810 × 540 2 cells and the nuclear energy release was boosted by factors of 3 × 10 3 , 1 × 10 4 , and 3 × 10 4 to study the dependence of the results on these parameters. Results. The boosted energy input results in convection at Mach numbers in the range of 10 −3 –10 −2 . Standard mixing-length theory predicts convective velocities of about 1.6 × 10 −4 if no boosting is applied. The simulations with AUSM + -up show a Kolmogorov-like inertial range in the kinetic energy spectrum that extends further toward smaller scales compared with its non low-Mach variant. The kinetic energy dissipation of the AUSM + -up solver already converges at a lower resolution compared to AUSM B + -up. The extracted entrainment rates at the boundaries of the convection zone are well represented by the bulk Richardson entrainment law and the corresponding fitting parameters are in agreement with published results for carbon shell burning. However, our study needs to be validated by simulations at higher resolution. Further, we find that a general increase in the entropy in the convection zone may significantly contribute to the measured entrainment of the top boundary. Conclusion. This study demonstrates the successful application of the AUSM + -up solver to a realistic astrophysical setup. Compressible simulations of convection in early phases at nominal stellar luminosity will benefit from its low-Mach capabilities. Similar to other studies, our extrapolated entrainment rate for the helium-burning shell would lead to an unrealistic growth of the convection zone if it is applied over the lifetime of the zone. Studies at nominal stellar luminosities and different phases of the same convection zone are needed to detect a possible evolution of the entrainment rate and the impact of radiation on convective boundary mixing. 

ABSTRACT Evolved Wolf–Rayet stars form a key aspect of massive star evolution, and their strong outflows determine their final fates. In this study, we calculate grids of stellar models for a wide range of initial masses at five metallicities (ranging from solar down to just 2 per cent solar). We compare a recent hydrodynamically consistent wind prescription with two earlier frequently used wind recipes in stellar evolution and population synthesis modelling, and we present the ranges of maximum final masses at core He-exhaustion for each wind prescription and metallicity Z. Our model grids reveal qualitative differences in mass-loss behaviour of the wind prescriptions in terms of ‘convergence’. Using the prescription from Nugis & Lamers the maximum stellar black hole is found to converge to a value of 20–30 M⊙, independent of host metallicity; however, when utilizing the new physically motivated prescription from Sander & Vink there is no convergence to a maximum black hole mass value. The final mass is simply larger for larger initial He-star mass, which implies that the upper black hole limit for He-stars below the pair-instability gap is set by prior evolution with mass loss, or the pair instability itself. Quantitatively, we find the critical Z for pair-instability (ZPI) to be as high as 50 per cent Z⊙, corresponding to the host metallicity of the Large Magellanic Cloud. Moreover, while the Nugis & Lamers prescription would not predict any black holes above the approx 130 M⊙ pair-instability limit, with Sander & Vink winds included, we demonstrate a potential channel for very massive helium stars to form such massive black holes at ∼2 per cent Z⊙ or below. 

ABSTRACT Massive stars are crucial to galactic chemical evolution for elements heavier than iron. Their contribution at early times in the evolution of the Universe, however, is unclear due to poorly constrained nuclear reaction rates. The competing 17O(α, γ)21Ne and 17O(α, n)20Ne reactions strongly impact weak s-process yields from rotating massive stars at low metallicities. Abundant 16O absorbs neutrons, removing flux from the s-process, and producing 17O. The 17O(α, n)20Ne reaction releases neutrons, allowing continued s-process nucleosynthesis, if the 17O(α, γ)21Ne reaction is sufficiently weak. While published rates are available, they are based on limited indirect experimental data for the relevant temperatures and, more importantly, no uncertainties are provided. The available nuclear physics has been evaluated, and combined with data from a new study of astrophysically relevant 21Ne states using the 20Ne(d, p)21Ne reaction. Constraints are placed on the ratio of the (α, n)/(α, γ) reaction rates with uncertainties on the rates provided for the first time. The new rates favour the (α, n) reaction and suggest that the weak s-process in rotating low-metallicity stars is likely to continue up to barium and, within the computed uncertainties, even to lead.